Giant resistance change across the phase transition in spin crossover molecules 
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The electronic origin of a large resistance change in nanoscale junctions incorporating spin 
crossover molecules is demonstrated theoretically by using a combination of density functional the- 
ory and the non-equilibrium Green's functions method for quantum transport. At the spin crossover 
phase transition there is a drastic change in the electronic gap between the frontier molecular or- 
bit als. As a consequence, when the molecule is incorporated in a two terminal device, the current 
increases by up to four orders of magnitude in response to the spin change. This is equivalent to a 
magnetoresistance effect in excess of 3,000 %. Since the typical phase transition critical tempera- 
ture for spin crossover compounds can be extended to well above room temperature, spin crossover 
molecules appear as the ideal candidate for implementing spin devices at the molecular level. 

PACS numbers: 



The last few years have witnessed several attempts at 
implementing spin devices at the single molecule level 
[T]-[3]. These include spin- valves, where organic molecules 
are used as transport medium for spin-polarized electrons 
injected from magnetic metals [IHB], and three termi- 
nal junctions incorporating magnetic molecules [7J. The 
second strategy appears particularly intriguing as the 
electrical response of single molecule magnets depends 
sensitively on both their spin and charge configuration 
[7J |H] . This sensitivity, combined with the possibility of 
manipulating the spin-state of a molecule by electrical 
means [9j [10] , makes spin-electronics based on magnetic 
molecules an attractive platform for both conventional 
and quantum logic. 

There are however two unfortunate drawbacks in con- 
structing spin-devices based on molecular magnets. The 
first concerns with the typically poor structural stability 
of such molecules away from solution. For instance al- 
most all the members of the Mni2 family react on metal- 
lic surfaces, so that their structural and electronic in- 
tegrity in uncertain once they are incorporated in a de- 
vice [TT] . The second problem has to do with the mag- 
netic anisotropy. In fact, although the anisotropy den- 
sity per atom may be very large, the overall magnetic 
anisotropy of the molecule is rather small, since only a 
handful of atoms contribute to the magnetic moment. 
This means that at room temperature the spin quan- 
tization axis of the molecule fluctuates at a frequency 
much higher than the transport measurement time and 
the molecule appears paramagnetic to a steady state 
transport experiment. Importantly, while the first prob- 
lem can be kept under control by engineering the end 
groups binding to the substrate [12] , the second one ap- 
pears much more tough to overcome. In fact, quantum 
tunneling of the magnetization is always present so that 
the instability in the magnetization direction may persist 
down to low temperatures. It is then not surprising that 
the only demonstration to date of spin- valve effect orig- 
inating from molecular magnets is limited to ultra-low 



temperatures [T3} fT4] . 

One possible way to overcome such a difficulty is to 
abandon the concept of spin- valve and to look at devices 
where the information is contained in the actual spin 
magnitude and not in its direction. A natural choice for 
such alternative device concept is that of using molecules 
displaying a magnetic bi-stability, i.e. molecules that can 
be found in two different spin states, both accessible with 
an external stimulus. This is the realm of spin-crossover 
compounds [15]. Such materials are formed by 3d tran- 
sition metal ions in an octahedral surrounding, which 
display a spin transition from a low spin state (LS), usu- 
ally the ground state, to high spin state (HS), usually a 
metastable one. Such a LS-HS transition can be triggered 
by temperature, pressure or light. It is accompanied by 
a modification of the geometrical structure, which alters 
the crystal field strength. The most extensively investi- 
gated spin crossover compounds contain Fe(II) and the 
transition is between a 1 Ai g LS state to a 5 T2 9 HS one. 

To date electron transport experiments in spin 
crossover molecules at the single molecular level remain 
scarce, mostly because the effort in depositing such 
molecules on surfaces has intensified only recently. How- 
ever, several deposition techniques are becoming avail- 
able (e.g. Ref. [16]) and a first demonstration of light- 
induced spin-crossover in vacuum-deposited thin films 
has been already reported [17]. Intriguingly the criti- 
cal temperature for spin-crossover in such thin films is 
similar to that of the bulk, giving hope that room tem- 
perature device operation may be soon achievable. Also 
intriguing is the possibility to tune the spin crossover 
transition with an electric field, as demonstrated recently 
for valence tautomeric compounds [T8] . 

Importantly there are two demonstrations of the inter- 
play between the spin-crossover transition and the elec- 
tron transport properties of a device. Prins et al. [19] 
measured a change in the current- volt age, I-V, curve of 
a Fe-(trz)3 (trz=triazole) complex deposited between Au 
nano-gap electrodes, which correlates well with the spin- 



crossover transition. Furthermore, the three-terminal 
single molecule experiments of Meded et al [20] showed 
low-energy features attributable to spin-crossover and 
controllable by gating. 

A crucial aspect of spin crossover compounds is that 
at the LS to HS phase transition there is a simultaneous 
electronic and structural change of the molecule. The 
electronic gap between the highest occupied molecular 
orbital (HOMO) and the lowest unoccupied molecular 
orbital (LUMO) varies accordingly to the different occu- 
pation of the 3d multiplet and at the same time there is 
an out-breathing relaxation of the octahedral cage con- 
nected to the weakening of the crystal field. As such, 
in transport experiments it is difficult to separate effects 
arising from the electronic structure from those originat- 
ing from the molecule geometry. In this letter we address 
such an issue with state of the art electron transport cal- 
culations based on density functional theory (DFT). In 
particular we show that the change in HOMO-LUMO 
gap at the phase transition may generate extremely large 
changes in the I-V characteristic. This is demonstrated 
for the prototypical Fe(II) SCC [FeL2] +2 , where L is a 
2, 2':6, 2 // -terpyridine group [see Fig. [TJa)]. Note that 
qualitatively similar results (not presented here) have 
been obtained for a second [FeL2] +2 compound where 
L = 2, 6— bis(pyrazol-l-yl), so that we expect our main 
conclusions to bare generality. 

Calculating the electronic properties of spin crossover 
molecules with DFT is a rather difficult task as the rela- 
tive energy separating the LS and HS state is very sensi- 
tive to the choice of exchange and correlation functional 
and to the atomic relaxation [21] (note that the molecule 
atomic coordinates are known only for single crystals but 
not in vacuum or on a surface). We then use the following 
strategy. We first perform structural relaxation for the 
molecule in vacuum by using the B3LYP functional [22] , 
one of the better performing for this problem. Then the 
results are compared with those obtained by using the 
generalized gradient approximation (GGA) and the local 
density approximation incorporating the self-interaction 
corrections (LDA-SIC) [23]. From this analysis we select 
the GGA functional to be used for the electron transport 
calculations, as this reproduces qualitatively the HOMO- 
LUMO gap obtained by B3LYP [24]. 

We start our analysis by discussing the electronic and 
geometric properties of the molecule in the gas phase. 
B3LYP calculations are performed with the NWCHEM 
package [25]. The basis set is 6-31++G** for H, C and N, 
while for Fe we use the Los Alamos double zeta with an 
effective core potential [26]. In Fig. [TJc) we present the 
total energy as a function of the reaction coordinates, 
X, interpolating between the LS (X = 0) and the HS 
(X = 1) geometry for the two different spin configura- 
tions. In the same figure we also show the dependence of 
the HOMO-LUMO gap, AE g , as a function of X. 

The figure clearly shows that the ground state is, as 





FIG. 1: (Color on line) Electronic properties of [FeL2] +2 , 
L = 2, 2':6, 2 // -terpyridine, along the reaction coordinates, 
X. In (a) and (b) we show respectively the 2,2 / :6, 2"- 
terpyridine group and the cell used for the transport cal- 
culations. Color code: C=orange, N=dark blue, H=white, 
S=green and Au=yellow. In panel (c) we present the total 
energy, E (full symbols to read on the left-hand side scale), 
and the HOMO-LUMO gap, AE g (open symbols to read on 
the right-hand side scale), for both LS and HS along X. Cal- 
culations are performed with the B3LYP functional. 



expected, low spin. As X increases there is a transi- 
tion to high spin for X^0.6 with the difference between 
the energy minima of the two spin configurations being 
about 200 meV. The HOMO-LUMO gap also changes as 
a function of X, increasing for both the spin states as one 
goes from the low to the high spin geometry (X gets big- 
ger), i.e. as the ligand field weakens. Most importantly 
at the relative energy minimum the gap, calculated with 
B3LYP, of the LS state is 3.87 eV (for X=0), while that 
of the HS is only 3.18 eV (for X=l). At the same two en- 
ergy minima LDA-SIC [27] gives us HOMO-LUMO gaps 
of 3.35 eV (LS) and 2.84 eV (HS), while GGA returns 
1.81 eV (LS) and 0.33 eV (HS). Although different ex- 
change and correlation functionals yield sensibly different 
gaps, in all cases this contracts significantly across the 
phase transition. Importantly both B3LYP and GGA 
give us spin down HOMO and LUMO levels in the HS 
state, while with LDA-SIC the HOMO is spin up and 
the LUMO is spin down. As such, we decide to use GGA 
instead of LDA-SIC for the transport calculations. 

We construct a two terminal device by attaching the 
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FIG. 2: (Color on line) Zero-bias transmission coefficient as 
a function energy, T(E; V = 0), for both the LS (red dashed 
line) and the HS (solid black line) configuration of the FeLi2- 
based junction and the molecular orbitals associated to their 
HOMO and LUMO like-orbital. Calculations are performed 
with the GGA functional. 



FeL 2 molecule to a Au(lll) surface via thiol group. In 
particular we consider bonding to the (111) hollow site 
and we optimize the Au-S distance for two different an- 
gles between the molecule and the normal to the surface, 
namely and 30°. As the transport properties depend 
weakly on such tilting angle the results reported here 
refer only to the 30° case, the bonding geometry pre- 
senting a lower energy. Transport calculations are per- 
formed with the Smeagol code [28 , 29 , which implements 
the non-equilibrium Greens function method within DFT 
[30] . The wave-function is expanded over a numerical 
atomic orbital basis set and the core electrons are de- 
scribed with norm-conseving pseudopotentials including 
non-linear core corrections. The basis set has double zeta 
quality for C, N, H and Fe, while only single zeta is em- 
ployed for Au. The simulation cell contains the molecule 
and 5 (lll)-oriented Au layers with 8x7 cross section. 
We use periodic boundary conditions in the direction or- 
thogonal to the transport with a uniform 2x2 fc-point 
sampling and an equivalent mesh cutoff of 400 Ry. The 
charge density is integrated over 64 energy points along 
the semi-circle, 64 energy points along the line in the 
complex plane and 64 poles are used for the Fermi dis- 
tribution (the electronic temperature is 25 meV). During 
the finite-bias calculations we integrate the Green's func- 
tion over the real axis on a 256 point mesh. 

Figure [2] displays the total zero-bias transmission coef- 
ficient, T(E; V = 0) = Y?* T<T { E ; V = 0), as a function 
of energy for both the LS and HS configurations of the 
junction (<j labels the spin, <j \). The most striking 
feature is a radical change of the position of the various 
resonant peaks as the molecule transits from LS to HS. 
This is a direct consequence of the re-arrangement of the 
Fe d shell density of state originating from the reduction 
of the ligand crystal field. In particular as the crystal 
field weakens in the HS state the energy width of <i-shell 
reduces and the peaks become more dense. This effect is 
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FIG. 3: (Color on line) I-V curve for a FeL2 molecule at- 
tached to gold electrodes in either the HS (black solid line) 
or the LS (red dashed line) state. The current is plotted on a 
logarithmic scale only for positive currents. The inset shows 
the spin- crossover magneto-resistance (SCMR) ratio, i?scMR. 
Calculations are performed with the GGA functional. 



particularly dramatic around the Fermi level, Ep, where 
the transport HOMO-LUMO gap shrinks. Recalling here 
the fact that in the linear response the conductance G is 
simply given by G = ^T(Ep] V = 0), we can conclude 
that at the phase transition there is a significant conduc- 
tance enhancement. One can quantify such an effect by 
defining the spin-crossover magneto-resistance (SCMR) 
ratio as i^scMR = (Chs — Gls)/Gls- At zero-bias we 
find i^scMR = 200 %, a variation which is certainly mea- 
surable. Note that across the phase transition there is no 
change in the molecule valence charge, so that the effect 
reported here occurs at an energy scale lower than that 
needed to observe Coulomb blockade. 

Some more details are provided by the insets of Fig. [2j 
where we show the frontier molecular orbitals associated 
to the different peaks in T(E; V = 0). In general around 
Ey the orbitals responsible for the transmission have all 
Fe-d character and mix little with the ligands. In the 
LS configuration the HOMO and LUMO are non-spin- 
polarized and have respectively d z i and d n character. In 
contrast the HS state has a HOMO of x 2 — y 2 symmetry, 
while the LUMO is c^-like, but both the HOMO and 
the LUMO are spin down. Thus the gap reduction is 
the result of the orbital re- arrangement of the d manifold 
associated to the spin crossover transition. In the LS 
state the gap is between the ti g triplet and e g doublet 
split by the octahedral crystal field. In the HS config- 
uration the system undergoes Jahn-Teller distortion in 
the singly occupied spin down triplet and thus the gap 
reduces. Note that such a distortion, and thus the spin 
crossover transition, may be obstructed if the molecule 
is too rigidly bound to the electrodes. As such for this 
effect to be detected one has to protect the spin active 
part of the molecule by appropriate chemical groups, so 
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FIG. 4: (Color on line) Transmission coefficients as a function 
of energy for different values of the gate potential V g . Here 
T(E; V = 0) is plotted as a color code on a logarithmic scale 
with red corresponding to low transmission (10 -11 ) and blue 
to high (10 _1 ). In the inset we show a schematic figure of the 
device investigated. The purple atoms represent the region in 
space where the gate is applied. The color code for the atoms 
is the same as in Fig. [TJb) . 



that the Fe-N cage can deform freely 

Then we look at the I-V characteristics, Fig. [3] 
Clearly, the drastic gap reduction in the HS state makes 
its relative current considerably larger than that of the 
LS. At a voltage of about 0.5 Volt the two differ by more 
than three orders of magnitude. The LS current is essen- 
tially tunneling-like over the entire bias window inves- 
tigated and it increases monotonically with bias. Such 
an increase is the result of the enlargement of the bias 
window, i.e. of the fact that a larger portion of the trans- 
mission spectrum contributes to the current. It is then 
the tail of the LUMO to contribute the most. The HS 
situation is different since voltages of the order of 0.4 Volt 
are enough for bringing both the HOMO and the LUMO 
transmission resonances within the bias window (note 
that orbital re-hybridization under bias drifts both the 
HOMO and LUMO close to each other hence decreas- 
ing the gap). As a consequence the current grows fast 
at low bias, when there is a transition from tunneling to 
resonant transport, and then saturates. Such a large dif- 
ference in the I-V's of the two spin states translates in 
a large i^scMR, which is presented as a function of bias 
in the inset of Fig. [3] Clearly .Rscmr as large as 3000% 
appear possible. 

Finally, we look at the effect of a potential gate on the 
transport properties of the molecule. This is applied by 
shifting the on-site energies of the molecular orbitals [31] . 
We note that for both states the transmission spectrum 
is characterized by a number of relatively sharp peaks 
spaced by at most 1.5 eV (most typically the spacing is 
~0.3 eV) and corresponding to states with substantial 



Fe d character. These are expected to shift with bias. 
As an example we consider the LS situation and in fig- 
ure gwe show T(E; V = 0) as a function of gate. Here 
T is plotted as a color code on a logarithmic scale with 
red corresponding to low transmission (10 -11 ) and blue 
to high (10 _1 ). Clearly the entire spectrum shifts with 
the gate potential until a state, either the HOMO or the 
LUMO, crosses the Fermi level (V = 0). At that point 
charge transfer takes place and the gate becomes less ef- 
ficient. Crucially when one of the two frontier molecular 
orbitals is brought to Ep the zero-bias conductance in- 
creases dramatically. Since the two spin states of the 
molecule have different frontier molecular orbitals, their 
position can be manipulated independently, i.e. they can 
be brought at Ep for different values of the gate. As such, 
an extremely large gate-modulated magnetoresistance is 
expected. 

In conclusion we have demonstrated that the changes 
in electronic structure associated to the spin crossover 
phase transition, namely the reduction of the HOMO- 
LUMO gap, produce large changes in the transmission of 
a two-terminal molecular junction. At low bias we have 
predicted a spin crossover magnetoresistances of the or- 
der of 200%, which can increase up to 3000% at finite 
bias. The effect is related to the different nature of the 
frontier molecular orbitals of the different spin configu- 
rations and can be tuned with a gate potential. Such 
an effect, whose experimental evidence has been already 
partially provided, can be exploited in the fabrication of 
spin-transistors at the molecular level. Intriguingly our 
results not only prove that spins can be manipulated in 
an organic molecule, but also introduce a mechanism for 
molecular switching which has electronic origin and it is 
not associated to a change in the molecule geometry [32] . 
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